Looking for an extra dimension with tomographic cosmic shear 



Yong-Seon Son^ 
Enrico Fermi Institute, University of Chicago, Chicago IL 60637 
(Dated: February 2, 2008) 

The cosmic acceleration was discovered in one of the brane-based models. We are interested in 
discriminating this model from the dark energy by tomographic cosmic shear. Growth factors are 
different in the two models when one adjusts parameters to get nearly identical H(z). The two 
models could be distinguished with independent determinations of both geometrical factors and 
the growth factor. We introduce new parameterizations to separate the influence of geometry and 
the influence of growth on cosmic shear maps. We find that future observations will be able to 
distinguish between both models. 
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I. INTRODUCTION 



The discovery of cosmic acceleration 0, 0, has led 
to theoretical efforts to understand the nature of the cos- 
mic acceleration 0, 0, IS an d phenomenological efforts 
to discriminate the diverse modelings [H, |tj 0, 0, 
IT3I Il4| . Most theoretical explanations for the acceler- 
ation preserve Einstein gravity and add a new smooth 
and dark component called Dark Energy (DE) 0, 
However, DE is not yet strongly supported theoretically 
nor phenomenologically. The similar cosmic acceleration 
was also discovered in one of brane-based models by 
C. Deffayet et.al. 0] (hereafter DGP model rep- 
resents this discovery of the cosmic acceleration in the 
brane-based model p.) 

In the DGP model, the cosmic acceleration is gener- 
ated by the modified gravity due to the presence of the 
extra dimension in DGP. The study of the extra dimen- 
sion physics influen cing the DGP model has been inves- 
tigated by 0, 0, Il9l |. Dvali, Gabadadze and Porrati 
designed a brane where ordinary matter is embedded in 
an infinite extra dimension and gravity is modified at 
larger scale @. The model (DGP) which generates the 
cosmic acceleration based upon |6|-type brane model was 
proposed in [HHItl . 

While the DE growth factor is suppressed by dark en- 
ergy component domination over matter component, the 
DGP growth factor is suppressed by the weakened grav- 
ity. We derive the DGP growth factor by using the set 
of continuity equations of ordinary matter confined to 
the DGP brane. We find that the DGP growth factor 
departs noticeably from the DE growth factor even with 
nearly identical expansion rate H(z). A slight difference 
is noticed from previous work based upon a different idea 

The constraints from cosmic shear maps on DE with 
a model de pen dent parametrization have been studied in 
ElEllllllllllliilll. The cosmological parameters 
related to DE are determined by a distinct evolution of a 



geometrical factor and a growth factor. We introduce a 
more general parametrization to determine directly, from 
cosmic shear data, the distance and the growth factor as 
a function of rcdshift. Both quantities can be precisely 
determined by the combination of the future CMB ex- 
periments and the future cosmic shear surveys. We in- 
vestigate how to probe the difference between the DGP 
growth factor and the DE growth factor with H(z) deter- 
mined through its geometrical influence on shear data. 



II. DGP MODEL 

We briefly review the DGP model and set up a formal- 
ism to draw our result. We consider a (4+l)-dimensional 
model in which no energy sources are present at an infi- 
nite bulk dimension except at a (3+l)-dimensional brane 
embedded in the bulk. The Einstcin-Hilbert action is 



S, 



(5) 



16tt 



(1) 



where the subscript (5) denotes that quantities are five- 
dimensional, and M(5) is the Planck mass in (4+1)- 
dimension. We add a (3+l)-dimensional brane with mat- 
ter fields and the induced metric to the action Eq. ^ The 
embedded three-brane action is [g 



5(4) = 



M 2 
16tt 



R 



(4) 



(2) 



where the subscript (4) denotes that quantities are four- 

(3+1)- 



dimensional, and M(a\ is the Planck mass 
dimension. 

We introduce a five dimensional metric with a line el- 
ement 



ds 2 



(5) 



= -n 2 (t, y)dt 2 + a 2 it, y)5 l jdx l dx 3 + dy 2 



(3) 



'Electronic address: ysong@cfcp.uchicago.edu 



where a flat metric is used for spatial dimensions and 
the coordinate y represents the extra dimension. Our 
universe is located at hypersurface y = where a lapse 
function is given by n(t, 0) = 1 and a spatial expansion 



2 



factor a(t, 0) is determined by the cosmic expansion of 
our universe [ill. 

In the DGP model, the gravitational force at short 
distances smaller than the crossover length r c scales as 
1 /r 2 and the gravitational force at large distances bigger 
than r c scales as 1/r 3 . The crossover length scale r c is 
given by 



(4) 



When the Hubble horizon is close to r c , the gravity 
is weakened and the cosmic expansion is accelerated. It 
is an alternative explanation of the cosmic acceleration 
with the absence of a dark energy component. We inves- 
tigate the phenomenological consequences of a DGP-type 
extra dimension in the following subsections, such as the 
geometrical factor and the growth factor. 



A. The geometrical factor of DGP 
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The Einstein equations from Eq. ^ and Eq. [21 are 



G 



(5)B 



R 



(5)B 



2 R 9(5)B 



^ rpA 



(5) 



The tensor T A B is composed of the energy momentum 
tensors of the bulk and the brane and the scalar curvature 
of the induced metric introduced in Eq. [3] We assume 
that the bulk is empty and the brane has only prcsureless 
matter components. Then the total energy p is 



3Af, 2 4) a 2 

P = Pm 1 <%)~2- 

sir ar 



(6) 



The hrst term of right side of Eq.[S]is the matter density 
confined to the brane and the second term comes from 
the scalar curvature of the brane 

By considering the first integral of Eq.[S]and the proper 
junction condition defining the derivative of a(t,y) in 
terms of y at crossing the surface y — 0, we get the 
following brane-FLRW equation 



FIG. 1: Probability distribution of fi m by fit the supernova 
data: top panel shows the fit of ACDM and the bottom panel 
shows the fit of DGP. 



where the constant tt c is fl c = (l/2r c Ho) 2 . The flat- 
ness constraint is imposed to relate Q c to fl m as Q c = 
[(1 — fl m )/2] 2 . Accordingly the crossover length scale 
is always greater than the present Hubble horizon to 
avoid the negative value of tt m or the empty universe 
(r c > H~ r ). There have already been attempts to fit 
the supernova data with DGP [23, El. W e use the 
'golden set' of supernova measurements |2Qj and com- 
pare the likelihood function of ACDM model with DGP. 
There is only a slight difference in the minimum x 2 be- 
tween ACDM and DGP, xLn = 178 for ACDM and 
Xmin — 179 for DGP. As far as the supernova experi- 
ment is concerned, both models are almost equally prob- 
able. ACDM model prefers high f2 m around 0.3 and DGP 
prefers low fl m around 0.2 (Fig. [TJ. 



tj2 H 8tt 



(7) 



where e denotes signature (+, — ). To be consistent with 
the observed cosmic acceleration, we choose the nega- 
tive sign for e. The negative sign implies de Sitter-like 
brane embedding in the extra dimension. The effect of 
the gravity leaking into the extra dimension is significant 
at H ~ l/r c but is negligible at H > l/r c pLIiq. 
From Eq. 0we get 



B. The growth factor of DGP 

We choose the Gaussian normal longitudinal (GNL) 
gauge to fix the degrees of gauge freedom in 5D gravity. 
In GNL gauge, the extra dimension is not perturbed and 
any transition of energy from the extra dimension to time 
coordinate or spatial coordinates is not allowed. Then we 
have the usual (3+I)-dimensional longitudinal gauge on 
the brane which leads to the perturbations in the time 
and spatial coordinates as 



H(a) = H 



fir 



i2 ) ; 



(8) 



500 = -l-2*(t,f,y = 0) J 

9ij = a 2 (t,y = 0)[l + 2<t>(t,x,y 



0)] k 



(9) 
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The 00 component of perturbed Einstein equation is 
given by [Hill 



2 r 

7^2 



fc 2 $ + 3d 



d - 



6 V M ( 3 5) 



P 6p, (10) 



where p is given by Eq. Eland 8p is given by 
8tt 2 



Sp = 5p n 



3M ( 2 4) a 2 



/c z $ + 3d 



d - d*^) 



(11) 



Neglecting the contribution of velocity perturbations, we 
write the Poisson equation by using Eq. Eq. ^3 an d 
Eq.Elas 



k 2 <£ = 



■in 



M, 



2-a 2 p„ 



1 



(4) 



1 



1 u rrn 



(12) 



2r c H 



where S m denotes the density contrast Sp m /p m . We do 
not consider e = 1 case which is ruled out by the current 
measure distance. 

The presence of Weyl component j32 does not deform 
the Poisson equation Eq. El The contribution of Weyl 
component is linearly added to the matter component 
with suppression factor of l/2r c H at early time. It could 
be important at later time when 1 /2r c H is close to unity. 
But in GN coordinate gauge, the trace of Weyl compo- 
nent vanishes and the fluid behaves like the radiation 
component. If Weyl fluid behaves like the radiation com- 
ponent then the clustering of that component is supressed 
and eventually it gaurantees the approximation of small- 
ness of Weyl contribution against the matter clustering. 

In the DGP model, matter does not flow into the ex- 
tra dimension, i.e. there is no momentum component 
along the extra dimension y. The perturbation length 
scales of interest in this paper range approximately from 
lO^ 4 !!^ 1 to lO -3 -?^ 1 . In this range, the gravity of den- 
sity fluctuations is not dissipated into the extra dimen- 
sion significantly since the physical scales of those modes 
are much greater than l/r c . Despite the presence of the 
extra dimension the set of coupled continuity equations 
of matter is not altered. In those intermediate scales, the 
coupled equations are 



drj 
drj 



+ tkVra = 0, 



aHv m — ik<&, 



(13) 



where we differentiate the matter perturbation S rn and 
the velocity perturbation v m in terms of the comoving 
time rj. 

Second order differential equation from Eq. 1131 can be 
written as 



By using Eq. [7] and Eq. ^1 the curvature perturbation $ 
is given by the background expansion rate H(z) as 



fc 2 $ = -a 2 H 2 - 
2 1 - 



l 

r c H 



(15) 



2r r H 



The curvature perturbation of DE is suppressed by the 
ratio of matter component to the gross energy in the 
universe, but the curvature perturbation of DGP is sup- 
pressed by the weakened gravity. The difference in causes 
of suppression leads to the distinct evolution of growth 
factors of DE and DGP. 



III. MODEL INDEPENDENT 
PARAMETRIZATION 

In this section we introduce a new parametrization to 
measure directly the geometrical factor and the growth 
factor by tomographic cosmic shear. First, we briefly 
discuss cosmic shear. 

The size and shape of galaxies is altered by gravita- 
tional lensing. The effect of lensing is described by the 
convergence k measuring the magnification or demagni- 
fication and the shear compone nts. 71 and 72, quantify- 
ing the distortion of shape [U 0| . The shear compo- 
nents can be inferred from measurement of galaxy ellip- 
ticities which are composed of an intrinsic ellipticity and 
a lensing-induced ellipticity. In the absence of correla- 
tions between the intrinsic ellipticities, the rms error in 
the measurement of each shear component is 



CT (7l) = cr(7 2 ) = 7rm 8 /V^ 



pix, 



(16) 



where 7 rms is the rms intrinsic shear of the galaxies and 
iVpjx is the number of galaxies in the pixel. 

We assume that photometrically-determined redshifts 
can be used to sort the source galaxies into eight redshift 
bins with Az = 0.4 from z = 0.0 to z = 3.2^23- Maps 
of the shear components can be decomposed into even 
parity E mode and odd parity B modes. The signal con- 
tribution to the covariance of shear E modes is correlated 
across bins, but otherwise diagonal. The E modes shear, 
7_e, power spectra are 

2 



drrWifu^Wifi^AKkx^), (17) 



where fi denotes the angular diameter distance to the 
median redshift of bin i. The geometrical factor appears 
in the window function W and the growth factor appears 
in dimensionless curvature power spectrum A| . We show 
how to measure directly those quantities in the following 
subsections. 



Parametrization of the angular diameter 
distance 



d 2 S d8 



drj 2 drj 



aH— - fc^$ = 0. 



(14) 



The lensing-induced ellipticities are weighted by the 
mean location of sources in the window function W(f , r) 
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TABLE I: The specifications of CMB experiments 



given by 

W(fi,r) = — (for r < fi, otherwise). (18) 

nr 

We directly parametrize the geometrical factor f, instead 
of determining f \ by the cosmological parameters depend- 
ing on the model choice between DGP and DE. 

B. Parametrization of the growth factor 

The dimensionless curvature power spectrum in Eg. 1171 

is 

A|(fc, r(z)) = ^-A g*(z) 2 T(k, z) 2 , (19) 

where Aq is the primordial amplitude, <?$(z) is growth 
function of <f> and T{k, z) is the transfer function. 
Here we are interested in the intermediate scales from 
lO -4 !^ 1 to lO -3 //,^ 1 where the perturbations of the 
dark energy do not grow significantly to influence on the 
matter perturbations in DE, and the gravity of pertur- 
bations does not dissipate into the extra dimension in 
DGP. In those scales, the transfer function with negligi- 
ble time dependence can be determined without the de- 
tailed knowledge of DE or DGP, and the departure from 
the linearity of density perturbations is not significant. 

We parametrize the growth factor by discretizing its 
evolution in each z bin as 

Fi = A ^, (20) 

where the growth factor Fi of matter perturbations is 
averaged inside each bin i and is normalized to Aq in the 
matter dominated era. The curvature growth factor g$> 
is related to the matter perturbation growth factor g by 
Eq.QJ 

IV. EXPERIMENTS AND RESULTS 

We sort out the cosmological parameters into three 
groups: primordial parameters, intermediate parameters 
and low redshift parameters. The primordial parameters 
come from a generic parametrization of structureless ini- 
tial conditions. The intermediate parameters contribute 



Experiment f sky n to t n to t/7o 
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0.5 


65 


1900 
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0.072 


100 


2920 



TABLE II: Weak lensing experimental parameters assumed. 
Units for the total source sky density ntot is 1/arcmin 2 and 
the per-component rms shear rms is evaluated at z — 1. 

to leave a signature on CMB acoustic peaks at last scat- 
tering surface. The low redshift parameters affect the 
density fluctuations at later times. The primordial pa- 
rameters consist of the primordial scalar amplitude Aq, 
the scalar spectral index n$ and the scalar running of the 
spectral index as- The intermediate parameters consist 
of the matter component u> m , the baryon component w^, 
the helium fraction yu e and the angular extent of sound 
horizon at the last scattering surface 8 s ■ The low redshift 
parameters consist of the neutrino mass m„, the reion- 
ization redshift z rc and all other parameters related to 
the cosmic acceleration. 

We consider the combination of CMB experiment with 
weak lensing survey in order to simultaneously constrain 
all cosmological parameters. The CMB angular power 
spectrum leads to tight constraints on the primordial pa- 
rameters and intermediate parameters. The power of 
constraints on all cosmological parameters with CMB 
alone is shown in |35| , and the constraints on the cosmo- 
logical parameters by the combination of CMB and weak 
lensing are extensively studied in |l3l l25j . The primor- 
dial parameters and intermediate parameters are tightly 
constrained with the future CMB experiments (see spec- 
ifications in TABLED). 

The primordial amplitude Aq can be precisely deter- 
mined by CMB alone [35J. The overall amplitude shift 
of CMB power spectra also results from the scattering 
of photons by ionized gas in the intergalactic medium. 
Though both cosmological parameters Aq and z rc are de- 
generate in most scales, the degeneracy between Aq and 
z rc is broken by the reionization bump of CMB polar- 
ization anisotropies in large scales. But the precision is 
limited by the unknown reionization history 36]. The 
primordial scalar amplitude Aq can be determined up to 
1.7 percent level of accuracy in the absence of lensing 
contribution. 

For the future weak lensing surveys, we assume uni- 
form coverage which gives a diagonal noise covariance 
matrix N, 

■^Imi. I'm' j ~ $ll'$mm'$ij (^-0 

n, 

where rii is the number of source galaxies per pixel in the 
given bin i. We use a ground-based survey G27t, and a 
space-based surveys S3000. The specifications for those 
surveys are shown in TABLE[n] The advantages of space 
are not fully exploited in this analysis which is only using 
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information on large-scale. The galaxy number distribu- 
tion for the ground-based survey G2tt are inferred from 
observations with the Subaru telescope ■ The analytic 
expression for this distribution is well-matched with 

dn/dz ~ z 13 exp [-(z/1.2) 1 ' 2 ] (z < 1) 
dn/dz ~ z ls exp [-(z/1.2) 1 - 2 ] (z > 1) (22) 



"(z)= 



We normalize the total number of galaxies n to t as 65 per 
sq.arcmin after reducing dn/dz a half in the 1.2 < z < 2.5 
due to inaccurate photometric redshift observation |38j . 
For the space-based survey, the distribution function is 



dn/dz 



1 exp(— z/1.5) 



(23) 



with reaching higher limiting magnitudes and having rel- 
atively accurate photometric redshifts |39|. 

The cosmological parameters considered in this pa- 
per are (A , ns, a s , w m , w b , Vh<s, &s, z xe , m v ) and (fi, ... 
,f$,Fi, ... , Fg). The Fisher matrix analysis is used to es- 
timate the constraints on all 25 cosmological parameters 
simultaneously. We show the constraints on fi and Fi in 
the following subsections. 



A. Constraints on the angular diameter distance 

The extensive window function defined by W(f, r) x 
T(k) is bell-shaped. As r is close to fi, the extensive 
window function decreases due to the numerator term 
of W in Eq. And as r is close to zero, the extensive 
window function also decreases since the transfer function 
decreases at high k = l/r. The variation of fj contributes 
to shift the peak point of the extensive window function. 
Larger median distance leads to more contribution of the 
smaller k modes to cosmic shear correlation functions. 

As is shown in Fig. |3 the angular diameter distance 
evolution is well-reconstructed by tomographic cosmic 
shear. The mean angular diameter distances fi are di- 
rectly measured in high accuracy around a few percent- 
age level. 



B. Constraints on the growth factor 

The variation of Fi has no impact on the shear correla- 
tions in bin j less than i. The response to the variations 
of Fi in bin j = 1 is a unique signature of a parameter F\ . 
Such a ladder structure helps to break degeneracy among 
the eight Fi parameters. The constraints on Fi at higher 
redshift bin are weakened since fewer shear correlation 
power spectra are influenced by Fi. 

Based upon the measurement of the primordial scalar 
amplitude Aq at the last scattering surface, cx(lnAo) ~ 
0.17, we proceed to construct the growth factor from 
higher redshift bins. The growth factor history is pre- 
cisely reconstructed up to z ~ 2 with the ground-based 
surveys and up to z ~ 1.2 with the space-based surveys 
as shown in Fig. [3J 
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FIG. 2: Constraints on the direct measurement of the angular 
diameter distance: The error boxes are the level of accuracy 
to determine fi in each bin. The solid curves represent the fi 
of Model DEr and the dash curves represent the fi of Model 
DEF. 
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FIG. 3: Constraints on the direct measurement of the growth 
factor: The error boxes are the level of accuracy to determine 
Fi as defined in Eq. HOI The solid curves represent the Fi 
of Model DEr and the dash curves represent the Fi of Model 
DEF. 
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V. CAN WE DISCRIMINATE DGP FROM DE? 

The direct measurements of the eight fj parameters 
and the eight Fj parameters allow us to find the pa- 
rameters that fit fj and Fj. In the DGP model, the 
fj and Fi are completely fixed by 9 other parameters 
(A , n s , a s , uj m , u b , y He ,9 s ,z rci m„) and the flatness con- 
straint. We assume that the fixed fi and Fj of DGP are 
fiducial values. In the DE model, the fj and Fi are still 
variable by the equation of state of dark energy w(a) even 
with all constraints given by 9 other parameters and the 
flatness. We parametrize w(a) asw+w a (l-a) |9j and 
vary (w,w a ) of DE to fit the fiducial fi and Fj. Model 
DEr denotes the best DE model to fit the fiducial fi and 
Model DEF denotes the best DE model to fit the fidu- 
cial Fi. In Fig. Hand Fig. El Model DEr and Model DEF 
provide the solid curves and the dash curves respectively. 




- 1 .0-0.9-0.8-0.7-0.6- 1 .0-0.9-0.8-0.7-0.6-0.5 



w w 



(w, w a ). However, any one of DE and DGP models is not 
excluded more than 90% confidence level with fitting fj 
alone. Also considering the generic aspect of DE, we can 
differently parametrize w(a) to fit better the fiducial fj. 
Thus it is not clear yet to discriminate DGP from DE by 
fitting f i alone unless DGP is excluded by its poor fit to 
n- 
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FIG. 5: Contour plots showing constraints on the DE param- 
eters with S3000, assuming the true model DGP. The contours 
in left panels cover DE parameter space to fit fi and the con- 
tours in right panels cover DE parameter space to fit Fi. The 
fiducial model is the DGP model and we vary DE parameters 
(w,w a ) to fit the fi and Fi. The inner curves represent 1-a 
confidence level and the outer curves represent 2-a confidence 
level. 



FIG. 4: Contour plots showing constraints on the DE param- 
eters with G27r, assuming the true model DGP. The contours 
in left panels cover DE parameter space to fit fi and the con- 
tours in right panels cover DE parameter space to fit Fi. The 
fiducial model is the DGP model and we vary DE parameters 
(w,w a ) to fit the fi and Fi. The inner curves represent 1-a 
confidence level and the outer curves represent 2-a confidence 
level. 

While the Fi of Model DEF is nearly identical with 
the fiducial Fj, the fj of Model DEr is distinct from the 
fiducial fj. The angular diameter distance to the last 
scattering surface at z ~ 1100 is almost fixed by the 
prior information of 8$ and u> m . When we fix the angu- 
lar diameter distance at z ~ 1100, we do not find any 
DE model which generates fj nearly identical with the 
fiducial fj by varying (w, w a ). By fitting fj, DGP is able 
to be discriminated from DE models parametrized with 



Next we consider fitting the f j and Fi simultaneously. 
Despite a good fit with the fiducial Fj, Model DEF is 
excluded by a poor fit with the fiducial fj. The distinct 
suppression in the DGP growth factor leads to the sep- 
aration in two regions of the DE parameter space to fit 
the fj and the Fj. In Fig. 0] and Fig. |3J we show the 
disagreement between two DE parameter spaces in de- 
tail. The ground-based cosmic shear survey G27r clearly 
discriminate DGP from DE around 95% confidence level 
(Fig. . And the space-based cosmic shear survey S3000 
is able to see the difference between DE and DGP around 
68% confidence level (Fig. [SJ . Although the space-based 
surveys in our analysis do somewhat worse because of the 
smaller sky coverage, we must remind the reader that we 
are focusing on science that can be done with informa- 
tion on the larger angular scales. Space-based observa- 
tions will have advantages with respect to ground-based 
observations at smaller angular scales. 

In conclusion, DGP is distinguishable from DE by a 
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cosmic shear survey. The DGP model is not yet clearly 
ruled out by any experiment. For instance, it is too early 
for us to rule out DGP by using a a at this moment. As 
shown in Fig. |3 the normalization influenced by DGP 
could be around 5% less. However considering the cur- 
rent constraints on the primordial amplitude Aq (more 
than 5% uncertainty), we are not able to discriminate 
DGP from DE by as alone with current datasets. We 
still have a chance to detect the extra dimension in the 
future. 

The DGP growth factor predicted by (2£j excluded the 
DGP model by more suppression than ours. Lue et. al. 
pof derived the gravitational perturbations of spherically 
symmetric clustered mater sources on the cosmological 
background. However the current solution is possibly 
missing the significant contribution of quadrature con- 
tribution (private communication with Roman Scocci- 



marro) and we are not able to compare both different 
methods at this moment. 

In case the future cosmic shear surveys prefer DGP 
to DE, then we will also be able to precisely determine 
the crossover length scale r c . The knowledge of r c will 
lead to the discovery of a fundamental energy scale and 
mark the beginning of experimental exploration of extra 
dimensions. 
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